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A general formulation of boundary conditions for semiconductor-metal contacts follows from a 
phenomenological procedure sketched here. The resulting boundary conditions, which incorporate 
only physically well-defined parameters, are used to study the classical unipolar drift-diffusion model 
for the Gunn effect. The analysis of its stationary solutions reveals the presence of bistability and 
hysteresis for a certain range of contact parameters. Several types of Gunn effect are predicted 
to occur in the model, when no stable stationary solution exists, depending on the value of the 
parameters of the injecting contact appearing in the boundary condition. In this way, the critical 
role played by contacts in the Gunn effect is clearly stablished. 



PACS: 05.45+b, 72.20.Ht, 85.30.Fg, 05.70.Ln 



I. INTRODUCTION 



The Gunn effect is an ubiquitous phenomenon in many 
semiconductor samples presenting negative differential 
resistance and subject to voltage bias conditions jiJ-Q. 
In a nutshell, the negative differential resistance makes 
it possible the existence of a variety of pulses and wave- 
fronts, which may be stabilized by the bias condition. 
Then a periodic shedding of waves by the injecting con- 
tact results in the periodic oscillation of the current 
through an external circuit, which constitutes the sig- 
nature of the Gunn effect. There is a vast literature 
on this topic, despite of which different basic questions 
concerning the Gunn effect remain poorly understood. 
Paramount among these, there are the questions about 
which are the correct boundary conditions and, given 
these, how to describe all the stages of the Gunn oscilla- 
tion. The lack of a precise formulation of the boundary 
conditions imposed by contacts on semiconductors and of 
a simple analytic treatment to analyze the Gunn oscilla- 
tions, has not allowed to clarify in a precise way the role 
played by contacts in the Gunn effect. It is worth noting 
that clarifying this point would open, for instance, the 
possibility of extracting information about the contacts 
from the analysis of the Gunn oscillations themselves, a 
subject of considerable interest for applied researchers. 



Recently progress has been made towards answering 
reasonably these two questions. On the one hand, ideas 
from Irreversible Thermodynamics || have been used 
to derive satisfactory boundary conditions for metal- 
semiconductor and other contacts in a general way |}]-^) . 
Previously the usual boundary conditions used in drift- 
diffusion semiconductor models were: (i) periodic, j^] 
(ii) charge neutrality, [Q (iii) fixed field, and (iv) 
control current-field characteristics of the contact |ll| ] 
plus phenomenological assumptions such as the "contact 
length" S. As boundary conditions (b.c.) for a semi- 
conductor presenting the Gunn effect, these conditions 
rank from clearly wrong (no current oscillation appears 
if the b.c. are periodic) to unsatisfactory because of their 
ad hoc character. Thus even when numerical simulations 
display the Gunn effect, the question of whether these 
results describe the real physical system where different 
contacts are present, is usually raised. In this paper, we 
shall present a simple derivation of appropriate b.c. for 
an ideal metal-semiconductor (MS) contact and use them 
to analyze the Gunn effect in Kroemer's model for bulk 
n-GaAs. Our description makes it clear which part of the 
derivation follows from general principles and which part 
includes the input from the physics of contacts. 

Concerning asymptotic descriptions of the Gunn effect 
which delve deeper than just numerical simulations of 
drift-diffusion models, some progress has been made re- 
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cently |12j-|l5[ . A detailed treatment of this topic can be 
found in Ref. @. 

The rest of the paper is as follows. In Section II we 
present our derivation of b.c. for ideal MS contacts and 
briefly discuss some other possibilities. Kroemer's model 
and its stationary solutions for these b.c. are analyzed in 
Section III. It is found that bistability between stationary 
solutions is possible for certain bias ranges depending on 
the values of certain dimensionless contact parameters 
io, ao, which are a combination of its effective density 
of states, barrier height, Richardson's constant, doping 
and temperature. Different types of Gunn effect, namely, 
charge monopoles (moving charge accumulation and de- 
pletion layers), and charge dipoles, (high and low field 
solitary waves), are predicted to appear depending on 
these contact parameters, when no stable stationary so- 
lution exists. In Section V we discuss our results, whereas 
Appendix A is devoted to technical matters related to the 
main text. 



II. BOUNDARY CONDITIONS 



The aim of this section is to present a systematic pro- 
cedure to derive b.c. at semiconductor contacts, estab- 
lished in previous works As a general rule, the 
method applies for non-degenerate semiconductors under 
moderate temperatures, that is, when thermionic emis- 
sion is the dominant transport process at the contact. 
Hence several contacts of interest, like ideal and non- 
ideal metal-semiconductor (MS) contacts or any type of 
heterojunction contact, can be modeled. Depending on 
the material parameters both limiting as well as ohmic 
contacts may then be described. It is worth noting that a 
precise modeling of this type of contacts may help to clar- 
ify the role played by other types of contacts used in semi- 
conductor systems, e.g. those in which thermionic-field 
or field emission processes dominate |l7| |l9j , for which a 
precise description, in the sense of the present paper, is 
not yet available. For the sake of clarity, the method will 
be presented along with its application to the case of an 
ideal MS contact. Other contacts have been considered 
in previous papers 



Let us consider an ideal MS contact. Due to the pres- 
ence of the contact, the magnitudes describing the phys- 
ical properties of the system, e.g electron density, elec- 
tric field, electron energy, . . . , may be discontinuous at 
that point. In addition, singular contributions localized 
at the contact itself, e.g. electron density at interface 
states (when they are present), may also occur. As a 
consequence, a given physical magnitude, d(x, £), can be 
decomposed as follows 

d(x,t) = d n (x,i)9(x) + d m (x,i)9(-x) +d s (t)S(x), (1) 



where d n , d m , d s refer to the values in the semicon- 
ductor (n), metal (m), and surface (s) parts , respectively 
(when no singular contribution is present, d s vanishes). 
Moreover, 9(x) is Heaviside's unit step function and 6(x) 
Dirac's delta function. They are introduced in order to 
represent the discontinuity across the contact and the sin- 
gular contributions, respectively. In writing Eq.(|l]), a one 
dimensional description of the system has been assumed, 
with the contact being located at x = and the metal 
[semiconductor] on its left [right] . By means of this type 
of decomposition, b.c. can be systematically derived. 

Our procedure consists of two steps: a) the identifica- 
tion of the relevant magnitudes describing the transport 
processes through the contact, and b) the derivation of 
precise laws describing such processes, which relate the 
relevant magnitudes at the contact and which constitute 
the desired b.c. For the first step, use will be made of a 
phenomenological formulation of transport through semi- 
conductor junctions Q while for the second of Shockley- 
Read-Hall (SRH) statistics |oLf|J. 



Let us consider a given magnitude, d(x, t), satisfying a 
standard balance equation of the form H 



dd(x,t) dJd(x,t) 



di 



di 



(Xd{x,t), 



(2) 



where J c i{x,t) and ad(x,t) refer to the current and net 
rate production associated to the magnitude d, respec- 
tively. It is not difficult to show that if similar balance 
equations were to be satisfied on each side of the junction 
and if surface fluxes only exist along the interface (what 
in a one dimensional description means Jd,s(t) = 0), then 
the following balance equation should be satisfied at the 
contact, [E2i 



dd s (t) 
di 



+ [Jd,„(0,i) - Jd,m(0,t)] = a d<s (t). 



(3) 



Now, we can proceed to calculate the net rate of en- 
tropy production at the contact, which will allow us to 
identify the relevant magnitudes describing the transport 
processes trough the contact. To begin with, we consider 
the balance equation for the total energy of the system. 
As this is a conserved quantity, we simply have 



de s (t) 
di 



+ [J e ,«(0,t)- J e , m (0,£)] =0. 



(4) 



As for an ideal MS contact no interface states are present, 
and hence no net charge or mass is accumulated at the 
contact, the total energy at the contact coincides with the 
surface internal energy, u s — e s p2fl . Hence, the balance 
of internal energy is described directly through Eq.(|^) or 
alternatively through 
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du s (t) 
dt 



+ [J„,n(0,i)- J«,m(Oj)] 



(5) 



with <r,, s (t) = [J u ,„(0,f) - J e , n (0,t)] - [J«,m(0,t) - 
>/e,m(0,i)]. In the previous expression, we have intro- 
duced explicitly the flux of internal energy (equivalent to 
the heat flux), which is in general different from the flux 
of total energy. Furthermore the Gibbs equation for an 
ideal contact is Tds s = du s (no interface states are 
present), where s s is the surface entropy and T the tem- 
perature. By assuming the contact to be in local equilib- 
rium, one then has Tds s (t) /dt = du s (t) / dt, which, after 
using Eq.(||), gives rise to the balance equation for the 
entropy 



ds s (t) 
dt 



+ [J s , n (0,i) - J s , m (0,t)] = C7., a (t), 



(6) 



with the entropy production given by 

<J s , s (t) = [J s , n (0,i) - -J e , n (0,f)] 

-[J s , m {0,t) - -J e , m (0,i}}. 



(7) 



A more explicit expression for tr s s is obtained once the 
bulk expressions for the fluxes are introduced on the r.h.s. 
of Eq.(Q). These expressions can be found elsewhere [SL||. 
One has, 

J s ,a(0,t) = ^J e ,a(0,t) - ^E F , a (0,i)J a (0,t); a = m,n 

(8) 



where Ef,o. refers to the electron Fermi level (or chemical 
potential) and J a to the electron number density current. 
Substituting into Eq.(m), we simply have 



c 5 , s (i) 



[E F}n (0, i)J n (0, t) - E F , m (0, t)J m (0, t)} 



[E F<n (0,i) - E Fm (0,t)}J n (0,t), 



(9) 



where in the second line use has been made of the con- 
tinuity of the electron number density current at an 
ideal MS contact (this continuity follows from the cor- 
responding balance equation for the electron number 
density by imposing that no carriers are accumulated 
(n s (t) = 0) nor created (<r n , s (t) — 0) at the contact). 
The final expression to be used in what follows, is ob- 
tained by introducing the electron quasi-Fermi levels, 
F a (Q,t) = E F , a (0,t) - eV a (0,t). Here V a (0,t) is the elec- 
tric potential (which is continuous through an abrupt 
junction) and e > is minus the charge of the electron. 
We then arrive at the desired expression, 

v s , s (t) = -^[^(0,t) - F m (0,t)}J n {0,t). (10) 



Eq. ( [To| ) shows directly that the relevant magnitudes de- 
scribing an ideal MS contact are the electron flux (elec- 
tron current density divided by e), J n (0,t), and the dis- 
continuity in the electron quasi-Fermi levels, (F n (0,t) — 
F m (0, t)), which plays the role of "thermodynamic force" , 
Q . Both flux and force vanish at equilibrium, and we as- 
sume (in accordance with the basic tenets of Irreversible 
Thermodynamics |^]) that there is a relation between 
them. When the fundamental relation between flux and 
force is specified, this relation is exactly the sought-after 
boundary condition at the contact. 

The relation between J„(0, t) and (F„(0, i) - F m (0, t)) 
should involve more information about the physics of 
the contact. First of all, let us notice that the entropy 
production ( |l0| ) is formally equivalent to the expression 
corresponding to generation-recombination processes |25| ] 
(or in general, to any activated process, such as unimolec- 
ular chemical reactions || or surface adsorption |m]|), 
provided J n (0,t) is identified with the net rate of the 
process. From this comparison we then conclude that the 
transport through an ideal MS contact may be described 
as an elementary kinetic process of the form: q m # q n , 
where q n and q m represent the carriers in the semicon- 
ductor and metal, respectively, with the net rate of the 
process being equal to J„(0,£). The kinetics of such a 
process can be described, for instance, by means of SRH 
statistics. As it is well known, this description relates 
the kinetic rate of the process, in our case J„(0, t), to the 
affinity, in our case the difference in quasi-Fermi levels, 
in agreement with our former general treatment. As it 
is shown in Appendix we obtain the following rela- 
tion between the net rate of the process, J„(0,i), and 
the jump of the quasi-Fermi level at the MS contact, 
(F n (0,i) - F m {0,t)): 



J„(0,t) = A e 



o /3(F„(0,t)-F m (0,t)) 



(11) 



E 



Here Ao is a positive constant [see ( |A3| ) in Appendix 
J], and e</>° = Ec(0,t) — F m (0,t) is the contact bar- 
rier height. Ec(0,t) = E^ — eV n (0,t), is the electron 
energy, with E^, the bottom of the semiconductor con- 
duction band and V n the electric potential at the semi- 
conductor surface. Moreover, (3 — (kT)^ 1 , where k is the 
Boltzmann constant. Eq.(|i"l|) is the sought b.c. for the 
ideal MS contact. It can be written more explicitly by 
using the expression h — Nce^^ Fn ^ Ec \ which holds for 
non-degenerate semiconductors with h being the semi- 
conductor electron number density and Nc the effective 
density of states. We then have, 



N c e 



n(0, t) 



(12) 
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In the non-degenerate case, Ao only depends on T (see 
Appendix This result ends the derivation of the b.c. 
for ideal MS contacts which we will use for the rest of 
this paper. It is worth emphasizing at this point that, 
as mentioned at the beginning of this section, the pro- 
cedure we have sketched here for the ideal MS contact, 
can be applied to several other types of contacts. These 
include non-ideal MS contacts with the presence of in- 
terface states Q and unipolar or bipolar heterojunction 
contacts with or without interface states 0. Moreover, 
by adding a few assumptions one can handle non-abrupt 
contacts ||. 

Note that for a MS contact located at x = L, that is, 
with the metal [semiconductor] on the right [left] hand 
side of the contact, the corresponding b.c. is 



J n (L,t) = (Nce-^i 



n(L,t) 



(13) 



We can compare our result, Eq.(|l2]) for the ideal MS 
contact with the corresponding one reported in |2f| [see 
also p. 261 of Ref. p6|]. Then we can identify Ai = 
AiT 2 /e, i — 0, L, where Ai is the Richardson constant 
for the semiconductor in contact with the metal located 
at i = 0,L. Theoretically, Ai, and hence Ai, would de- 
pend only on the given semiconductor but not on the 
metal 26 . However, in practice Ai is taken as a phe- 



nomenological parameter and it can not only depend on 
the metal but also on the preparation procedures (27| . On 
the other hand, basic energetic arguments lead immedi- 
ately to the following rule for the contact barrier height 
^6),g8|: 4>b = 0m — X- Here (f> % M is the work function of 
the metal in the MS contact located at x = i and \ is 
the semiconductor electron affinity. For covalent semi- 
conductors, the validity of this rule has been put under 
question for the last five decades p7| , ^8[ . However, re- 
cently it has been shown that even for this type of semi- 
conductors, if accurately growth materials are used, good 
agreement is obtained with this simple rule p9j , ^0| . No- 
tice, that when such materials are not used, as very often 
happens, the contact formed turns out to be non-ideal. 
This is so because it is difficult to avoid that very thin 
insulating layers and/or interface states may be present 
at the contact [^8 19 1. Hence, a non-ideal description 
of the contact should be used, with for instance, a bias 
dependent relation for the barrier height which would 
include the effects of the insulating layer and of the in- 
terface states. Such contacts have been described in de- 
tail in Ref. |7]||. In particular, it has been shown that 
for this non-ideal contacts there is no simple description 
using directly equations such as ([ll]) or ([l2]), for con- 
tact non-stationary effects induced by interface sates (not 
present for ideal MS contacts) introduce additional terms 
not present in Eq.([TT|). These terms could be of impor- 
tance when describing naturally non-stationary phenom- 
ena such as the Gunn effect with non-ideal contacts, and 



will be considered in future works. 

Eqs. (|l|) and (|l3|) can be rewritten in terms of the 
electric field by using the Poisson equation to eliminate 
the electron density from them: 



9E . ~. e / „ T «, fl j< , J n (i,t) 



N C 



(14) 



Here the upper [lower] sign holds for i = [i = L]; e 
and no are the bulk semiconductor permitivity, and its 
doping, respectively. 

We have thus shown how our procedure allows us to de- 
rive explicit and precise expressions for the b.c. imposed 
by a given contact. A first important consequence of this 
method can be drawn directly from Eq. (14), which is 



simply a relation between the normal derivative of the 
electric field and the current density at the contact. Ex- 
amining our derivation shows that this result is simply 
a consequence of the use of kinetic models to describe 
the exchange of carriers through the contact. Hence, 
one should expect that b.c. derived in this way, will 
result in relationships between the normal derivative of 
the electric field and the current density at the contact. 
It is easily seen that (if diffusion effects can be neglected) 
our b.c. can be transformed into a Kroemer's type con- 
tact current-field control characteristics jll| (see next sec- 
tion). However, unlike previous models following Kroe- 
mer's approach ||,[I|,[l2],|l] , our control characteristics is 
the result of a physically precise derivation, and therefore 
only parameters which are physically well-defined appear 
in it. In particular, to use our control characteristics we 
do not have to invoke ad hoc assumptions involving new 
parameters such as the contact length ||,|l|. 

To facilitate the analysis in the rest of the paper, it 
is convenient to rewrite Eq. ( |l4|) in dimensionless units. 
This greatly reduces the number of relevant parame- 
ters. Our dimensionless electric field £ , electron densi- 
ties n, current densities j{x,t), time t and position x are 
measured in units of Eq, Uq, ehoHoEo, e/(e^o«-o) an d 
£E /(eno), respectively |52|. In these equations, E and 
Ho are an electric field and the zero-field electron mobil- 
ity typical of the processes occurring in the bulk of the 
semiconductor (see below). Then Eq. (113) becomes 



d£ 

— (i,t) = ai[-ii±j(i,t)], 



(15) 



where we have defined 
HoEqNc 



i 



N, 



no 



(16) 
(17) 
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As mentioned before, here the upper [lower] sign refers 
to i = [i = L]. It is worth noting that on is always a 
positive quantity (because \ is) while it does not have 
a definite sign: it depends basically on the value of the 
barrier height, (j>\, and on the doping value, tiq. 

It should be noted that there are important restric- 
tions on the possible values of the contact current density 
which are due to the fact that in Eqs. (|T^) and ( |l3| ) the 
electron density n, is a positive quantity: 



< h(i,t) = N c e 



Nc 



or in dimensionless units 

±j(i,t) <aj l +U=jt at . 



(18) 



+ U = if at is always a 



[Eqs. (|l|) and (0) imply that a, 
positive quantity, equal to the maximum current density 
which the contact can provide] . These restrictions on the 
current are reminiscent of the rectifying properties of MS 
contacts. In practice, they only impose a real limitation 
for the case of true rectifying contacts (when one of the 
jsat j g sma ii) Otherwise, i.e., for large values of jf at , an 
ohmic contact is obtained which does not impose a real 
limitation on the current. 

In order to analyze the influence of the derived b.c. 
on the Gunn instability, we will assume a sample formed 
by a certain semiconductor (able to display the Gunn ef- 
fect) and by two MS contacts implemented on it. The 
resulting b.c. are 



d£ 

— (0,t) = a \j(0,t)-i } 
d£ 



to^ 



-oti\iL +j(L,t)}. 



(19) 
(20) 



As discussed above, for a given semiconductor the values 
of the contact parameters may vary somewhat depending 
on the metal used in the contact and on the preparation 
procedures. For instance cti may vary two orders of mag- 
nitude, from about 0.3 to 33.4, if we use the experimen- 
tal values of Richardson's constant for GaAs reported in 
p7| . Similarly, the values of ii may also span two or- 
ders of magnitude, from about 0.03 to 4.01, due to the 
variation of oti, if we fix the barrier height, 4>b ~ 0.2 V 
(corresponding to Al |^9| ) , and the donor density is 10 14 
cm" 3 . Thus there is a rather wide range of parameter 
values for the contacts, corresponding to a large variety 
of situations which will be described in this paper. 

Lastly, the applied bias V, defined as eV{t) = 
F m (L,t) — F m (0,t), can be expressed as follows: V(t) — 
£(x,t)dx + (0° — <f)L). In the previous expressions, 
V and 4>\ are in units of eE^/efiQ and F m in units of 
eEl/fiQ. 



Very frequently the analysis of the Gunn effect under 
dc voltage bias is carried out by using the opposite sign 
for the electric field: E = —£, Then the b.c. become 



dE 

— (0 > t) = ao[*o-i(0 > t)] I 
dE 

— (L,t) = a L [i L +j(L,t)} 



(21) 
(22) 



With these definitions, the dc voltage is V = — j Q Edx — 
[c\>\ — 0°). Instead of working with the voltage V, it is 
more convenient to use the average electric field on the 
semiconductor sample, 4> = L~ 1 J Q L E(x, t)dx, which is 

equal to cf) = j;[—V + (0" — In what follows, neg- 

ative voltages, V < 0, will be considered such as > 0. 
With these conventions, the carriers go from the cath- 
ode (injecting contact) at x = to the anode (receiving 
contact) at x = L. 



III. KROEMER'S MODEL AND ITS 
STATIONARY STATES 



A. Kroemer's model 



The unipolar drift-diffusion model for the Gunn effect 
proposed by Kroemer i is generally accepted to pro- 
vide a rather complete description of the main features 
of this effect. Yet it is simple enough to allow very de- 
tailed asymptotic analysis; other important models such 
as Biittiker and Thomas's [Q incorporate more detailed 
physics but their study is technically more demanding. 
In the dimensionless units described above, Kroemer's 
model is 



dE (dE 

i 

L 



. d 2 E 
dx 2 



J. 



E(x, t) dx = (f>. 



(23) 
(24) 



Eq. ( |23|) is Ampere's law which says that the sum of dis- 
placement current and drift-diffusion current is equal to 
the total current density J(t) . It can be obtained by dif- 
ferentiating the Poisson equation, dE/dx = n — 1, with 
respect to time, substituting the charge continuity equa- 
tion dn/dt+dj(x, t)/dx = [the electron current density 
is of the drift-diffusion type: j(x,t) = nv(E) — 8dn/dx], 
and then integrating the result with respect to x. The 
electron velocity is assumed to be N-shaped and for spe- 
cific calculations we shall use [03| 



v{E) = E 



1 + BE 4 
1 + E* 



(25) 
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(it has a maximum vm > at Em > followed by a 
minimum < v m < vm at E rn > Em) and the electron 
difusivity S to be constant. The results using other curves 
having the same shape are similar. If v(E) does not reach 
a minimum but saturates instead as E — > oo, not all the 
monopole and dipole waves which we have found occur. 
Thus we have chosen the velocity curve that yields the 
richest dynamical behavior. The behavior of Kroemer's 
model with saturating velocity will be commented upon 
in the discussion. The dc bias <f> is the average electric 
field on the semiconductor sample. Eqs. (p3[)-(p4[) need to 
be solved with an appropriate initial field profile E(x, 0) 
and subject to the following mixed boundary conditions 
resulting from substituting j(x,t) = J(t) — dE/dt [from 
(1)] into ©-©: 

dE f OE \ 

— (0, t) = a U - J(t) + ^(0, t) , (26) 

^(L,t) = a L (i L + J(t)-^(L,t^ . (27) 



(In the case 5 — 0(1) the procedure is slightly more com- 
plicated and we shall omit the corresponding details; see 
p2"t . In this section we shall analyze the stationary states 
of Kroemer's model in n- GaAs'lJH under dc voltage 
bias with the new boundary conditions (pl|)-(p2]). Our 
work is based upon previous asymptotic and numerical 
studies of this and related models |l2|-|l5| . 

In this asymptotic limit, any stationary solution can be 
described as composed of outer and inner solutions: the 
outer bulk solution is a piecewise constant field profile 
valid everywhere except for two narrow boundary lay- 
ers located at the contacts and, for particular values of 
the current density, a narrow transition layer somewhere 
in the middle of the sample (see Fig. |l] and explanation 
below). First of all, if we ignore inner solutions, the sta- 
tionary state solves the equations 

v(E)-J = 0(e), 

£ = </> + 0(e). (31) 



In what follows, i, will be assumed to be positive be- 
cause the physically interesting phenomena (including 
the usual Gunn effect mediated by high field domains) 
are observed for these values of ii, as will be seen in the 
following sections. 

For typical n-GaAs data, S <C 1 and L » 1 Q. In 
this limit, we shall find approximate solutions to the ini- 
tial boundary value problem (^3|)-(^7|) for E(x,t) and 
J(t). Strictly speaking, the simple asymptotic descrip- 
tion that follows holds in the limit L — > oo, even when 
S = 0(1) |ljj. Assuming 5 <C 1 just simplifies the de- 
scription of the traveling waves of electric field in the 
semiconductor through the use of characteristic equa- 
tions and shock waves |l2] , |35| , |36|| . 

To take advantage of this limit, we will use the follow- 
ing rescaled time and length, 

1 t 

= T > V 



L 



x 
L 



Then Eqs. (§f)-(§§) become 



E (y, s ) dy = (f>- 



(28) 

(29) 
(30) 



except for particular values of J which will be speci- 
fied below. These equations result from retaining only 
order-one terms in Eqs.(|29|) and d30|), and assuming 
E(y) =const. Then for those values of <j> such that the 
outer solution ( [Si] ) is compatible with the boundary con- 
ditions, we have J = v(<j>) + 0(e). Let us denote by 
E\(J) < E 2 (J) < E 3 (J) the three zeros of v(E) - J 
\E 2 (J) is unstable]. Then the outer (bulk) field profile 
will be E(y) = Ei, i = 1,3, depending on the value of 
the bias <j>. 

At y = and y = 1 there are boundary layers, which 
we will call injecting and receiving layers, respectively. 
E(y; J), the field at the injecting boundary layer of width 
O(e) at y — 0, obeys [we ignore the O(S) diffusive term]: 



dE J 
ay v(E) 
0E 

e— (0; J) = a (i - J), 
dy 



(32) 
(33) 



The analysis of Eqs.(|32|) and (|33| ) is more easily carried 
out if we express the derivative b.c. Eq.(|33j) in terms of 
a b.c. for the electric field at the contact, E(0; J) . We 
can obtain E^O; J) from ( |33] ) by using ( |32| ) to eliminate 
8E(0; J)/dy. The result is that E = E(0; J) solves 



B. Stationary states and its stability 



ME) 
ME) 



J, where 
(1 +a i )v(E) 
l+a v(E) 



(34) 
(35) 



Before describing the Gunn effect in the present model, 
it is convenient to discuss how to construct the station- 
ary solutions of the model in the limit e<Cl and S <C 1. 



Notice that the contact curve j c (E) has the same ex- 
trema as the velocity curve v(E) and saturates for high 
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electric fields to the value jQ at , defined in Eq.([18|). Kroe- 
mer's contact characteristic for shallow-barrier metal- 
semiconductor contacts presented in Q , corresponds to 
a particular case of our model in which the electrons in 
the metal are assumed to be in equilibrium with those 
of the semiconductor near the contact. For this case, 
one would take ao — > with \otoio\ < oo, so that j c {E) 
would be then proportional to v(E). In contradistinction 
with Kroemer's contact characteristic, the general curve 
j c (E) may intersect the bulk velocity curve v(E). The 
main difference between these two cases is that a Gunn 
effect mediated by charge dipole waves is seen only if 
j c (E) intersects the bulk velocity curve v(E). If (|34j) has 
a solution, Eq. (|32|) indicates that E(y; J) approaches one 
of the solutions of ( |3ll) as we leave the boundary layer. 
The boundary layer at the receiving contact y — 1 is a 
much narrower diffusive boundary layer of width O(eS). 
The field there is l3ll 



i-y 

eS 



dE _ r E ' 
dv Je 



v(E) dE, 



V = 



e l v(E) dE 



1 f Ei 

where «l (il + J) = / v(E) 

6 Je,. 



dE, 



(36) 
(37) 

(38) 



</-</> characteristics are flat, with J = Jf for certain con- 
stant values of the current. The corresponding outer 
field profile is step-like with E = E 2 (Jf) = E c2 (Jf) 
if < y < AY and E = Ei(J f ) with i = 1,3 if 
AY" < y < 1 [see Fig. [lj. AY" is chosen so as to sat- 
isfy the bias condition <j) ps E 2 {Jf)AY+(l - AY)Ei(Jf). 
The flat part of the J-<j> characteristics corresponds to 
a bias range E\(Jf) < <f) < E 3 (Jf) (see the details be- 
low). Finally, when J is near the value jg a *, the field 
at the injecting contact is very large, the contact re- 
gion is almost depleted of electrons and its extension, 
y = e[Ecs(J) — Ei{J)] (i = 1,3), may be compara- 
ble to the length of the sample. Assuming the exten- 
sion of the depiction layer at the injecting contact is 
less than the sample length, the corresponding bias is 
4> « \e[E c3 {J) - E t { J)] 2 + Ei(J) (i = 1, 3). The charac- 
teristics tends to saturate at j'q '. 

Following this general scheme, different possibilities 
may be distinguished according to the relative values 
of j cm < jcM < jo at with respect to v m < v M - Here, 
jck = (1 + a i )v k /(l + a v k ) and v k , with fc = m,M, 
refer to the minimum (fc = m) and maximum (fc = M) 
of the contact, j c (E), and velocity, v(E), curves, respec- 
tively . We now discuss the different cases which appear 
for our velocity curve. 



(i = 1, 3) whenever ( |38[ ) has a solution El. 

The idea now is to fix J and to discuss for which values 
of J the above construction yields a stationary solution. 
Additionally, its stability will be considered. Clearly we 
may distinguish different cases according to the values of 
the contact parameters (ii, oti), i = 0, L. In what fol- 
lows we shall assume for the sake of simplicity that the 
boundary layer equations at the receiving contact, Eqs. 
( [37| ) and (j3g), always have a solution, and hence only the 
parameters of the cathode, (io, ap), need to be consid- 
ered. 

The general situation encountered when constructing 
the stationary solutions is the following. For each value 
of J, there are one or three values of the contact elec- 
tric field at x = 0, which are solutions of Eq. ([54]). We 
shall denote these field values by E C i(J), with E C \{J) < 
E C 2(J) < E c3 (J). Then the field profile in the inject- 



ing boundary layer is a monotonic solution of Eq. (32), 
which joins E c i(J) (i = 1, 2, 3) to one of the solutions of 
J = v(E) (outer solution). Furthermore, the outer solu- 
tion may be a constant field profile given by E(y) = Ei (J) 
(I = 1,3) which extends to the end of the sample, where 
a narrow receiving boundary layer exists [see Fig. |IJ . For 
such an electric field profile, the bias is <fi w Ei{J). The 
corresponding J-<fi characteristics satisfies J ~ v(<fi). By 
this construction, we identify the portions of the J-<f> 
characteristics which follow the first or the third branch 
of v(E) (see the details below). Other portions of the 



1. jcm < jcM < jo^ < Vm < VM 



In this case, see Fig. || , we have for < J < J c m a 
class of solutions joining E C ±(J) and E\{J), with voltage 
4> w Ei(J). For < 4> < Ei(j cM ) the curve J-<j> then 
follows the first branch of v(E). Furthermore, a second 
class of solutions joining E c j,{J) and E±(J), will exist for 
jcm < J < jo at - In this case for J not near j™', the 
voltage is given by <f> rs E\{J), and for J near j™' it 
is by ps ±e[E c3 (J) - E^J)} 2 + E X (J). Then, in the 
characteristics the third branch starting at cj> ps Ei(j cm ) 
follows the first branch of v(E) at the beginning, till it 
tends to saturate to Jq"' for larger voltages (see Fig. || ). 
Joining these two classes of solutions, there exists a third 
class for j cm < J < j c u which joins E c2 (J) to Ei(J), 
with (f> pa E\{J). These solutions are unstable, and they 
give rise to the second branch in the characteristics, Fig. 
||, which also tend to follow the first branch of v(E). 
Note that for voltages Ei(j c m) < <j> < E\{j c M), the two 
classes of stable stationary solutions coexist [see inset at 
the bottom of Fig. |2j. Hysteresis between them is then 
possible. 
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2. jcm <V m < jcM < jo"* < V M 

In this case, see Fig. || , the description is very sim- 
ilar to the previous case, except on what concerns the 
third branch of the J-tfi characteristics. Now, this branch 
is composed of two types of solutions: i) for j cm < 
J < v m there is a class of solutions joining E c3 (J) 
and Ei(J). Most of the time one has 4> ~ E\(J), ex- 
cept for J near v m that the solution is step-like with 
<j) w E m AY + (1 - AY)E 1 (v m ). We expect that these 
solutions become unstable on a bias range which is a 
subinterval of Em < 4> < E m |37|,|8). Then a Gunn ef- 
fect mediated by moving charge monopoles (which are 
charge depletion layers) might appear (see companion 
paper [^6|). For v m < J < j$ at , there is a class of so- 
lution joining E c3 (J) and E 3 (J), with (j> ss E 3 (J) for 
J not near j s at and (j> w §e[£ c3 (J) - E 3 (J)] 2 + E 3 (J) 
for J near Jq°* . Then, the third branch of the J-(f> 
curve starts following the first branch of the v(E) curve 
for Ei(j cm ) < 4> < E\{v m ), then it presents a flat re- 
gion for E\{v m ) < 4> < E m with J = v m , correspond- 
ing to the step-like solutions, and finally a region for 
E m < (ft which starts following the third branch of the 
v(E) curve and tends to saturate to jQ at for larger volt- 
ages. Note the presence again of bistability for voltages 
Ei{jcm) < (j> < Ei(j C M) [inset at the bottom of Fig. ||] 
and hence the possibility of hysteresis. 

A similar situation to the one depicted above would 
appear for j cm < j cM < v m < j^ at < v M - 



3. Vm <C Jcm < jcM < jo < V M 



moving charge dipoles) will be present for these values of 
the bias (see the companion paper |]l6| ). 

A similar situation appears for v m < j cm < j c m < 
v M < io at - 



4- v m < j cm < vm < jcM < jo at 



In this case (see Fig. ||), the third branch of the J-cf> 
characteristics is described as in the previous case. The 
first branch is composed of two types of solutions: one 
class, for < J < %, joining E cl {J) and Ei(J), and 
the other, for vm < J < j C M, joining E C \(J) and E^(J). 
For the first type of solutions, one has 4> ~ E\{J) except 
for J w v M that = E M AY + (1 - AY)E 3 (v M ). These 
step-like solutions are expected to become unstable in a 
subinterval of E M < <j) < E m |^7| , p8| . Then a Gunn ef- 
fect mediated by moving charge monopoles (which are 
charge accumulation layers) might appear (see the com- 
panion paper |l6|]). Thus, the first branch starts follow- 
ing the first branch of the v(E) curve for < (f> < Em, 
then it presents a flat portion for Em < <f> < E 3 {vm), 
with J — vm , and ends following the third branch of the 
v{E) for E 3 (vm) < <j) < E 3 {j c m)- The second (unstable) 
branch of the J-(f> curve is formed by a class of solutions 
that starts at E C 2(J) and ends at E 3 (J), with <f> ~ E 3 (J). 
Then, this branch follows the third branch of the v(E) 
curve, for E 3 (j cm ) < <\> < E 3 (j C M)- Note that for this 
range of bias, two stationary stable solution coexists [see 
inset at the bottom of Fig. ||] and hysteresis may appear. 



The main difference for this case with respect to the 
previous ones, relies on the second and third branches 
(see Fig. ^). Now, the third branch of the J-cf) curve 
involves one single class of solutions joining E c3 (J) and 
E 3 (J), with voltages </> « E 3 (J) for J not near j'o at , and 
(j) « \e[E c3 (J) - E 3 (J)} 2 + E 3 (J) for J near f Q at . The 
second (unstable) branch is formed of two classes of so- 
lutions: one class for i < J < j C M starting at E c2 (J) 
and ending at Ei(J), and another class for j cm < J < ig, 
starting at E c2 (J) and ending at E 3 (J). For J « i , these 
solutions are step-like with the voltage given through 
« E c2 (i )AY + (1 - AY)Ei(J), with i = 1, 3 depend- 
ing on the class of solutions considered. The (unstable) 
branch in the J-<j> curve then starts following the first 
branch of the v(E) curve for Ei{j c m) < <t> < Ei(i ), 
then it presents a flat portion for Ei(i ) < cj) < E 3 (i ) 
with J Ri i , and it ends following the third branch of 
the v(E) curve for E 3 (j cm ) < <p < E 3 (i ). Note that 
for voltages Ei(j c m) < <j> < E 3 (j cm ) there is no stable 
stationary solution [see inset at the bottom of Fig. [|. 
Thus we expect that the usual Gunn effect (mediated by 



C. Phase diagram 



By collecting the information obtained in the previ- 
ous subsections, the phase diagram describing the differ- 
ent behaviours of the system can be sketched, in terms 
of the injecting contact parameters, io and uq, Fig. 
^. Stable, non-oscillatory stationary solutions are ex- 
pected for values of io and ao such as jg at < v m , where 
j^ at = 1 + i . Otherwise, some kind of oscillatory so- 
lution should be found. Charge accumulation monopoles 
appear for j C M > vm (or equivalently for io > vm), 
charge dipoles for v m < j cm < j C M < "m, that is, for 
v m < io < vm, and charge depletion monopoles for 
jo at > v m with j cm < v m (i a < v m ). For complete- 
ness, also the separation between low field and high field 
dipoles, discussed in the companion paper [^6), has been 
depicted. It is worth noting that this rich phenomenology 
of oscillatory states appears just by changing the value 
of the contact parameters. This fact should be taken 
into account in analyzing the Gunn effect in real sys- 
tems, where, as mentioned before, a wide range of values 
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for the contact parameters, depending on the metal used 
and preparation procedures, may appear. 



IV. DISCUSSION 



We have presented a general formulation for the deriva- 
tion of the boundary conditions imposed by metal- 
semiconductor contacts on semiconductor systems. Ac- 
cording to this general formulation, the appropriate 
boundary conditions for ideal metal-semiconductor con- 
tacts are linear relations between the normal derivative of 
the electric field at the contacts and the electron current 
there. For the classical unipolar drift-diffusion Kroemer's 
model of the Gunn effect, these boundary conditions are 
of mixed type. In this paper, we have investigated how 
the boundary conditions for ideal metal-semiconductor 
contacts affect the stationary solutions of the Kroemer 
model, and their stability. Depending on the values of 
the contact parameters, bistability and hysteresis may 
appear. Moreover, for some range of parameters no sta- 
ble stationary solution is expected to occur. In those 
parameter ranges we expect to find the Gunn effect. Nu- 
merical simulations show that different types of Gunn 
effect appear, mediated by a variety of waves: (i) charge 
monopole accumulation wavefronts, (ii) monopole deple- 
tion wavefronts, or (iii) charge dipole waves (high and low 
electric field domains). Why these types of Gunn effect 
appear in the simulations will be explained by the asymp- 
totic theory of the companion paper It suffices to 
say that without this theory we would have missed signifi- 
cant possible instabilities. For example (ii) seems to have 
been missed by earlier workers, in spite of past extensive 
simulations of Kroemer's model ||. With our bound- 
ary conditions, the previously described types of Gunn 
effect are found in the following ranges of dimensionless 
critical contact currents: (i) jcM > vm (ii) jem < v m 
and jo°* > v m, and (iii) v m < j cm < j C M < vm- Here 
jcM = jc{E M ), jcm = jc(E m ), jo at are the critical cur- 
rents, and v m and % are the minimum and maximum 
values of the electron drift velocity v(E), E > 0. When 
we want to characterize experimental samples display- 
ing the Gunn effect, it is important to bear in mind the 
great influence of the contact parameters on the type of 
wave mediating the Gunn effect. A wide range of values 
for these contact parameters may be obtained depending 
on the type of metal used or on the contact preparation 
procedure followed. 

N-shaped velocity curves occur naturally in recently 
observed self-sustained oscillations in weakly-coupled n- 
doped GaAs/AlAs superlattices (see Ref. |59| for the 
most complete data so far). In these superlattices there 
is strong indirect evidence for a Gunn effect mediated 
by charge accumulation monopoles through photocurrent 
and photoluminescence measurements [flOf. It is hard to 



say at this point which other possibilities of those found 
in our analysis might be realizable in these systems. An 
important issue to be decided is the form of the boundary 
conditions. Our analysis needs to be modified in order 
to be extended to these systems, as quantum tunneling 
plays an essential role in the injection of carriers through 
contact regions. 

The most used velocity curves v(E) for the classical 
Gunn effect in bulk n-GaAs lack the third branch after 
v m . The reason is that avalanche breakdown appears 
at electric fields smaller than E m . The avalanche field 
is smaller for the longer samples needed to observe the 
Gunn effect and this precludes reaching the high fields on 
the third branch of v(E). Then low field dipole domains 
and charge depletion monopoles are not observed in the 
usual bulk samples or in strongly coupled superlattices 
with wide minibands, which are analogous to them p4]. 
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APPENDIX A: DERIVATION OF EQ. (0) BY 
MEANS OF SRH STATISTICS 



Let us consider the elementary kinetic process: q m ^ 
q n describing the charge transport through the junction. 
By assuming the validity of the SRH statistics to describe 
this process, the following expression for its kinetic rate, 
J„, can be obtained p3[ : 

Jn = J dE n J dE m D m (E m )D n (E n ) 

X[fm(E m )(l - fn(E n ))7mn(E m , E n ) 

-/«(£«)(! - f m {E m )) lnm {E m ,E n )], (Al) 

where D a (E a ), a — n,m, is the density of states of sys- 
tem a, f a (E a ) its occupation function, given through 
the FD distribution f a (E a ) = (1 + e 0(E a -E Fa )yl ( with 
Ef.u the corresponding Fermi level. ~f mn (E m , E n ) [resp. 
lnm{E m , E n )\ is the probability per unit time for the 
transition between states of energy E m and E n [resp. E n 
and E m ] . At equilibrium, we must have J n = and F m = 
F n , with F a = Ep, a — eV a being the corresponding quasi- 
Fermi levels. This implies Ep, n — Ep^ m = e (V n — V m ) = 
(using the assumed continuity of the electric potential at 







the contact). These equations follow from Eq. ( Al) if the 
latter is supplemented with the following detailed balance 
relation: 



-7 (E E 1=7 (E E ) P P( E ™- E i) 

lmn\- LJ 'm t -L^n J — ]nra\ J - i ra^ )^ 



(A2) 



A term (3 (V n — V m ) has to be added to the argument of 
the exponential in (A2) if V n ^ V m ; see (8) for a more 
general case. We now substitute Eq. (A2) into Eq. (Al) 
and use the equations 



E C = E° C - eV n , 



e<jP b =E C - F„ 



(Eq is bottom of the semiconductor conduction band and 
e§\ is the height of the contact barrier). After straight- 
forward manipulations, we derive (|ll]), in which the tran- 
sition coefficient Ao is: 

A = y" dE n J dE m D n (E n )D m (E m )(l - f n (E n )) 

x(l - f m {E m ))lnrn(E mi E n )e^ E o- E -l (A3) 

When the semiconductor is non-degenerate, we may ap- 
proximate 1 — f n {E n ) w 1, whereas for a metal we may 
approximate fm(E m ) by its equilibrium value. Then, for 
this case, Aq is a function of T only. 
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FIG. 1. Stationary electric field profiles, showing the piece- 
wise character of the solutions. The dashed line corresponds 
to a step- like stationary solution. Narrow boundary layers are 
present at the contacts. 



1.0 



0.8 



_ 0.6 



0.4 



0.2 



0.0 





0.11 



0.07 

0.07 0.11 0.15 



0.0 



0.2 



0.5 



0.8 



1.0 



FIG. 2. Stationary current-voltage characteristics, for 
L = 500, i = 0.048 and a = 9, for which 
jcm < jcM < jo ai < v m < dm, showing bistability for bi- 
ases Ei(j c m) < cf> < Ei(j c M). The dashed line corresponds 
to the unstable solutions with E(0;J) — E C 2(J). For com- 
parison the v{(j>) curve is also plotted (dotted-dashed line). 
Insets: At the bottom, blowup of the bistable region. On 
top, contact characteristics, jc(E), (dashed line) and velocity, 
v(E), (continuous line) curves for this case. 




FIG. 3. Stationary current-voltage characteristics, for 
L = 500, io = 0.135 and ao = 8, for which 
jcm < v m < jcM < jo at < "m, showing bistability for biases 
Ei(jcm) < 4> < Ei(j C M)- The dashed line corresponds to the 
unstable solutions with E(0;J) = E C 2(J). The flat portion 
of the curve corresponds to J = v m . A Gunn effect mediated 
by moving depletion charge monopoles is expected on a bias 
range which is a subinterval of Em < 4><E m - For comparison 
the v{4>) curve is also plotted (dotted-dashed line). Insets: 
At the bottom, blowup of the bistable region. On top, con- 
tact characteristics, j c {E), (dashed- line) and velocity, v(E), 
(continuous line) curves for this case. 




FIG. 4. Stationary current-voltage characteristic for 
L = 500, i = 0.27 and a = 4, for which 
Vm < jcm < jcM < Jo"' < "m, showing the unstable sta- 
tionary solutions (dashed line) with E(0; J) — E C 2(J). The 
flat portion of the curve corresponds to J = i . Note that no 
stable stationary solution exists for -Ei(jcm) < 4> < Es(j cm )- 
Then a Gunn effect mediated by moving charge dipoles is 
expected. For comparison the v(<p) curve is also plotted 
(dotted-dashed line). Insets: At the bottom, blowup of 
the unstable region. On top, contact characteristics, j c (E), 
(dashed- line) and velocity, v(E), (continuous line) curves for 
this case. 
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FIG. 5. Stationary current-voltage characteristic for 
L = 500, in = 0.87 and a = 0.5, for which 
*\ Jem < »m < jcM < jo at ) showing bistability for biases 
E3(j cm ,) < 4> < Es(jcM)- The dashed line corresponds to the 
unstable stationary solutions with E(0; J) = E C 2(J). The flat 
portion of the curve corresponds to J = v m . A Gunn effect 
mediated by moving depletion charge monopoles is expected 
on a bias range which is a subintcrval of Em < <f> < E m . 
For comparison the v((f>) is also plotted (dotted-dashed line). 
Insets: At the bottom, blowup of the bistable region. On 
top, contact characteristics, j c (E), (dashed-line) and veloc- 
ity, v(E), (continuous line) curves for this case. 
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FIG. 6. Phase diagram, showing the diversity of Gunn os- 
cillations that may appear depending on the values of the in- 
jecting contact parameters, in and an. The different separatri- 
ces correspond to jo at = v m (continuous line), j cm — v m (dot- 
ted line) and j C M = vm (dotted-dashed line). Also depicted 
the separation between low and high field charge dipoles, dis- 
cussed in the companion paper (dashed line) . 
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